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A magnetic helicity integral is proposed which can be applied to domains which are not magnet- 
ically closed, i.e. have a non-vanishing normal component of the magnetic field on the boundary. 
In contrast to the relative helicity integral, which was previously suggested for magnetically open 
domains, it does not rely on a reference field and thus avoids all problems related to the choice 
of a particular reference field. Instead it uses a gauge condition on the vector potential, which 
corresponds to a particular topologically unique closure of the magnetic field in the external space. 
The integral has additional elegant properties and is easy to compute numerically in practice. For 
magnetically closed domains it reduces to the classical helicity integral. 



INTRODUCTION 

Magnetic helicity is an important quantity in describ- 
ing the structure and evolution of magnetic fields in many 
fields of physics, in particular in plasma physics and 
astrophysics. It was introduced to Plasma Physics by 
H.K. MofFatt in and was originally defined as an inte- 
gral over a magnetically closed volume, i.e. a volume for 
which the normal component of the magnetic field (B) 
vanishes on the boundary dV: 

H{B):= [ A-BdV] Bn^B-n\ov = 0. (1) 
Jv 

Here A is the vector potential, V x A = B, of the mag- 
netic field. The integral measures - roughly speaking - 
the Gaussian linkage of magnetic flux within V. More 
precisely, it is the asymptotic linking number of pairs of 
field lines averaged over the volume |2| ■ It is an important 
property of this integral that we can derive an equation 
of continuity for the helicity density, which uses only the 
homogenous Maxwell's equations (here E = —dtA — V4> 
is the electric field): 

c't(A-B) + V-((/)B-t-E X A) = -2 E-B . (2) 

It can be shown that the integral is a topological invari- 
ant, i.e. it does not change under a deformation of the 
field within V , as given for instance by the motion of a 
magnetic field embedded in an ideal plasma, satisfying 
(v is the plasma velocity) 

E + V X B = . (3) 

Under such a condition, (jSJ becomes 

at(A-B) + V-((0- v-A)B + vA-B) =0 , (4) 

so that integrating over a volume with v-n = on the 
boundary (or more generally a comoving volume) results 



in 

4 / A-BdV = I 9t(A-B) + V-(v A-B) dl/ 

dt Jv Jv 

= - / (0- v-A)B-n(ia 0. (5) 

JdV 

Moreover, the total helicity is often an approximate in- 
variant for non-ideal plasmas Q , and is therefore a valu- 
able tool in determining the evolution of many technical 
and natural plasmas. One of the earliest results was the 
prediction of the relaxed state of a Reversed-Field Pinch 
Q, but there are many more applications (see for an 
overview) . 

However, the boundary condition _B„ = on the inte- 
gral which is necessary to ensure gauge invariance, 
means that it can not be applied to cases where the mag- 
netic field crosses the boundary. Typical examples are 
the vacuum vessels of laboratory plasmas where an ex- 
ternal magnetic field crosses the boundaries, or the at- 
mospheres of stars or planets, where the studied volume 
is usually bounded by the surface of the body, through 
which the magnetic field emerges. 

In such cases it was previously necessary to resort to 
the calculation of the relative helicity, i.e. the helicity was 
calculated with respect to a reference field Bj-ef satisfying 
the same boundary conditions. One can prove 0,0 that 
for an arbitrary closure of the magnetic field outside V , 
denoted by Boxt, the relative helicity 

if(B|Bref) = if(B + Bext)-i?(Brcf +Be.t) (6) 

= y"(A + Aref)-(B-Brcf) dT^ , (7) 

is actually independent of the external closure of the field. 
The reference field is in most cases choosen to be a po- 
tential field (see e.g. |1|) since a potential field is easy 
to compute and physically distinguished as the lowest 
energy state compatible with the boundary conditions. 
The introduction of a reference field, however, not only 
complicates the calculation of magnetic helicity, but also 
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complicates its already difficult interpretation. For in- 
stance, the question arises as to whether a change of 
relative helicity in a volume has a physical meaning, or 
whether it is only due to our particular choice of reference 
field. 

In this contribution it is proposed to replace the refer- 
ence field by a more general boundary condition on the 
vector potential and it is shown that this leads to a well 
defined quantity. 

DEFINITION 

The volume V C Si^ we refer to is assumed to be 
simply connected and without cavities, i.e. it has vanish- 
ing first and second Betti numbers. The first condition 
ensures that the vector potential is unique up to a gra- 
dient of a function, while the second implies that a vec- 
tor potential always exist (alternatively we can require 
J B ■ nda = over the boundary of any cavity) . More 
general volumes such as a solid torus can be considered 
as well, if additional constraints arc included. For a solid 
torus, for example, we have to impose J A-dl = around 
the hole of the torus. 

Definition: The universal magnetic helicity of a mag- 
netic field in a simply connected volume V C JR^ is de- 
fined as 



Hv{B) 



= A^A = on the boundary, which is a closed mani- 
fold, and therefore X\qv = constant. This means that A 
is constant on the boundary, but it will in general vary 
inside V. However, a gauge transformation A — > A + VA 
with X\dv = const, conserves the helicity integral: 

iJy(B) ^ Hv{B)+ [ VX-B-nda 
Jv 

/ VA-B = / XB-nda^X B-nda = 0(ll) 

Jv JdV JdV 

Existence: Starting with an arbitrary vector potential 
for B in y, with Vt-At = p, we note that it is possible 
to find a function A(ui, U2) on the boundary with A-^A = 
p. Here ui and U2 are boundary coordinates as defined 
above. This function can be extended to a function on 
all of V, for instance, by letting it smoothly fall off to 
zero within an e- neighbour hood of the boundary A(x) :— 
A(wi,U2)/(w3) with 



exp(-M§/(e2 - uj)) U3 < e 







elsewhere 



(12) 



A-B dV wither -At = on dV. (8) 



Thus A = A — VA has the desired property Vt-At ~ 0. 
Note that this proof also shows how to satisfy the bound- 
ary condition in practice. In particular there exist stan- 
dard numerical routines to solve A^A = p on arbitrary 
boundaries. 



Here Vt • At = is the divergence of the tangential 
component of A on the boundary dV, i.e. the divergence 
is taken with respect to the boundary coordinates only. 
For an explicite calculation we can choose an orthogonal 
curvilinear coordinate system (ui, U2, U3), defined locally 
on the boundary such that the unit vector U3 coincides 
with n . Then {ui,U2} span the tangent space of the 
boundary and Ay is represented as Ay = A^Ui -\- A^U2 
Denoting the scaling factors hi = ||^|j the divergence 
reads 



1 



d 



hih2 \dui 



h.A' 



d 
du2 



hiA' 



(9) 



In order for the quantity HyiB) to be well defined, we 
have to prove firstly that it is not gauge dependent, and 
secondly that the boundary condition can be satisfied for 
any field. 

Gauge invariance: First note that the curl of At on 
the boundary, that is in a two-dimensional surface is a 
scalar. Using the coordinate system from above it reads 
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Vt X At = -r^ [ 77-/12^2 



d 



4-hiA' 
0U2 



Ba . (10) 



Together with the boundary condition Vt ■ At = 0, it 
uniquely determines At, since a gauge A A + VA con- 
sistent with the boundary condition requires Vt • (VA)t 



PROPERTIES 

Showing that the universal helicity is well defined is 
not enough to justify its name. It must also reduce to 
the total helicity for the case of a vanishing normal 
component of B on the boundary. Since ((HJ includes the 
case of vanishing i?„ this is obviously the case. 

In addition, we can prove that it is a topological in- 
variant for any deformation of the magnetic field inside 
V, i.e. for any deformation, which leaves the boundary 
unaffected: v ~ 0\dv- Evaluating (jS)) on the boundary 
gives Et = —dtAr — {y<p)T = 0. Since At is deter- 
mined solely by B„, which is constant in time, we obtain 
(9f At = and hence </> = const, on the boundary. Thus 
(ISJ vanishes, now due to </> = const, and v = 0\dv instead 
of due to Bn ~ 0. 

Another important property of the universal helicity 
integral is its additivity with respect to magnetic fields. 
The rule is the same as for the total helicity: 

i/y(B'^ + B'') =i/y(B") + i/y(B'') + 2iJy(B^B''), (13) 

where 

iJy(B",B'')= I A°--B.^dV^ I A^-B" dV (14) 
Jv Jv 
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is the mutual or cross helicity integral. The equivalence 
of the two integrals in 1)14(1 can be shown by using the 
condition Vt-Aj^ = 0, which implies a representation 



At = V^E' X n. 



(15) 



In order to proof this note that Ay has a representation 
as At = slt x n with another vector field tangential 
to the boundary. Comparing (O with (|10|l we sec: 

= Vt ■ At- = Vt X ar ■ (16) 

Thus the field is a gradient. Analogous we find 

-B„ = Vt X At = -Vt • slt = -At* (17) 

An example of At and the potential 'I' on the boundary 
is shown in Fig. ^ 




FIG. 1: An example of two adjacent regions of positive and 
negative B-n = — At*!' on the boundary. The grey lines are 
contours of the potential ^' defined by 1151 . They coincide 
with field lines of At. In black field lines of Vvl* are shown. 

The representation H15|) can now be used to prove that 
the difference between the two integrals in 1] 14(1 vanishes: 

[ V-(A" X A'')dF = / (A"xA'')-nda (18) 

Jv JdV 



dV 



dV 



fV*" X V*'')-n da 



V X f^V*")-!! da = 0. 



Here the last integral vanishes due to the application of 
Stokes theorem to a surface without boundary. 

Furthermore, the universal helicity is additive with re- 
spect to complementary volumes. Complementary here 
describes that the volumes and are adjacent to 
one another, and that they satisfy = Bf^ on their 
common boundary, and have = B^^ = on all other 
boundaries. Note that this still assumes that the volumes 
are simply connected and have no holes. Thus, the total 
volume y° U has vanishing normal magnetic field on 
its boundary, and we can calculate its (classical) total 
helicity 




FIG. 2: A flux surfaces of the field B'' derived from ^ for 
the same boundary condition as used in Fig. 



The proof relies on the fact that the total helicity on the 
left hand side is gauge invariant, so that we can choose 
a gauge for the vector potential such that Vt • At — 
holds both on its boundary and on the interface between 
V"- and V''. Then the total vector potential can be split 
in two parts with A° = A|ya (A° = on V'^) and A** 
analogously defined. This implies (|19|l . 



INTERPRETATION AND EXAMPLE 

In the following we prove that the boundary condition 
Vt • At = corresponds to the existence of a unique 
field in the exterior domain V'' such that A'' ■B'' = 0. 
Thus the univeral helicity integral ffya(B'') equals the 
total helicity integral of completed by B** as Eq. l(T^ 
shows. In order to see this we use a coordinate system 
as introduced for the representation Q and define A** in 
a layer of thickness e from the boundary with /(us) as 
defined in fl^ : 



A** = At(ui,U2)./(m3) 



(20) 



B" = Br,iu,,u■2)f{u^)u' + ^^\/^{u^,U2)i21) 

One easily checks that this satisfies A'' • B'' = 0. Note 
that the field B** is tangent to surfaces spanned by U3 
and V^f. Such a surface is shown in Fig.|21 In the limit 
of e ^ the field B** becomes a singular " surface field" , 



B" 



S„(ui,U2)i?(-U3)u'-'5(«3)V*(Mi,U2) . (22) 



i7, 



(B" + B^) - i?yo(B") + Hy^iB'). 



(19) 



Here S{x) is the Dirac delta function and H{x) the Heavi- 
side function. That is, B** is a field of finite magnetic flux 
in the surface, which diverts the normal component Bn 
in a field along V^'. This field is uniquely determind by 
the boundary condition. 

Example. For a non-trivial example, consider a field 
consisting of two untwisted flux tubes with the same mag- 
netic flux $ as shown in Fig. Ol^a). There are no closed 
field lines in the volume under consideration, so there is 
no way of calculating the helicity of this configuration 
with the classical helicity integral. The configuration of 
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the singular boundary field is shown in Fig. |3fb). An 
explicit calculation yields 77y(B) = — <i>^. This can be 
understood as the helicity of the closed field shown in 
Fig-E^a), which shows a closure of the field with a config- 
uration topologically equivalent to the singular boundary 
field. The closed tube has a uniform twist of — 27r and 
thus a total hehcity of i?v(B) = 

Another way of calculating TJy (B) is to make use of the 
symmetry of the configuration and use Eq. H19|l . Fig. 01 
shows one example. The magnetic field is closed with 
an identical copy in V'. This configuration has a total 
helicity of —2$^ due to the linkage of the two flux tubes 
(see e.g. 0). Hence -2$^ = iJy,(B)+i?v(B)= 2Hv{B) 
and therefore i?v(B) = — $2. 




a) b) 

FIG. 3: a) The example under consideration and b) the cor- 
responding boundary with contours of (gray) and field lines 
of V* (black). 




SUMMARY 

In this letter it was shown how the total helicity inte- 
gral can naturally be generalized to allow for magnetic 
fields which are not closed within the domain, i.e. which 
have a non-vanishing normal component on the bound- 
ary. The construction does not require an explicit ref- 
erence field as the relative helicity integral does, which 
was previously used in this situation. Instead we have 
a gauge condition for A on the boundary which corre- 
sponds to closing the domain with a topologically unique 
field. This field is an external complement with zero he- 
licity density to the field in the given domain. The new 
integral has all desirable properties, i.e. it is gauge invari- 
ant, topologically invariant, and it reduces to the total 
helicity whenever the latter is well defined. Moreover, 
it shows the proper additivity with respect to fields and 
complementary volumes. This facilitates not only many 
calculations of helicity, but also its interpretation. 
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FIG. 4: a) The field shown in Fig. El is closed by two loops 
which are topological equivalent to the singular boundary 
field, b) A closure of the field with Hv'{B) = HviB). 



